Magnetic impurity resonance states and symmetry of the superconducting order 

parameter in iron-based superconductors 
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We investigate the effect of magnetic impurities on the local quasiparticle density of states (LDOS) 
in iron-based superconductors. Employing the two-orbital model where 3d electron and hole con- 
duction bands are hybridizing with the localized /-orbital of the impurity spin, we investigate how 
various symmetries of the superconducting gap and its nodal structure influence the quasiparticle 
excitations and impurity bound states. We show that the bound states behave qualitatively different 
for each symmetry. Most importantly we find that the impurity-induced bound states can be used 
to identify the nodal structure of the extended s-wave symmetry (S ± ) that is actively discussed in 
ferropnictides. 

PACS numbers: 74.20.Rp, 74.25.Ha, 74.70.Tx, 74.20.-z 



I. INTRODUCTION 

The problem of magnetic impurities in a su- 
perconductor has been extensively discussed in the 
literature 1 ^' 3 -^' 5 -. The magnetic impurity and its moment 
can interact with the conduction electrons of the metal 
in the normal or superconducting state. In the former 
case this leads to the Kondo effect and a resonance state 
at the Fermi level. In the latter case it is well known that 
a single magnetic impurity doped into a superconductor 
produces a localized bound state within the quasiparticle 
excitation gapi. The spectrum is sensitive to the sym- 
metry of the order parameter and is therefore a powerful 
tool to probe the pairing symmetry. 

The discovery of new Fe-based superconductors 6 - with 
distinct multi-orbital band structure 7 -!^ have opened a 
new horizon to high temperature superconductivity. One 
of the most significant questions for these materials is the 
symmetry of the superconducting gap and the underly- 
ing Cooper-pairing mechanism. The latter is believed 
to arise due to purely electronic mechanism and a vari- 
ety of models have been investigated with various weak- 
coupling approaches within random phase approximation 
RP A 10 : 11 ' 12 and renormalization group technique s! 13 ' 14 . 
It was concluded that the fully gapped extended s-wave 
state with the 7r-shift of the gap between electron and 
hole Fermi surface sheets is the most natural outcome 
of these theories. It is believed to be driven by the in- 
tcrband spin fluctuations at the antiferromagnetic wave 
vector (-7T, 7r) in folded Brillouin zone and it also competes 
with the spin density wave instability at the same wave 
vector which leads to the columnar or striped AF state 
for low doping. 

However, despite intensive experimental efforts, the 
pairing symmetry of this new class of superconduct- 
ing materials is not completely settled. Some ex- 
perimental groups have reported the fully gapped 
behavio r 15 : 16 : 17 : 18 : 19 , but some measurements, in partic- 



ular NMR relaxation and penetration depth suggest ex- 
istence of gap node a 20 ' 21 ! 22 . From the theoretical side 
it has also been realize d 11 ' 23 that the superconducting 
gap structure may be non-universal in ferropnictides due 
to the large intraband Coulomb repulsion. Its inclusion 
may force the superconducting gap to develop a node 
which crosses one of the Fermi surfaces. At the same time 
the symmetry of the superconducting gap will, however, 
still remain extended s-wave though higher harmonics 
are acquired. Moreover, in some scenario a 24 ' 25 the super- 
conducting gap even changes from the extended s-wave 
towards either d x 2_ y 2- or (i^-wave symmetries depend- 
ing on the slight variation of parameters. Recently it 
has been found that isoelectronic substitution of As by 
P in BaFe 2 (Asi_ 2 ,P K ) 2 changes the gap structure in Fe 
pnictide compounds from nodeless to nodali 6 -. It seems 
that whereas electron- and hole doping leads to a fully 
gapped state, isoelectronic doping (equivalent to chem- 
ical pressure) leads to the presence of line nodes in the 
gap function. Therefore it is desirable to investigate the 
magnetic impurity effect in Fe pnictide compounds for 
different candidates of gap symmetries. The resulting 
characteristics of the LDOS which is sensitive to the 
nodal structure may provide a clue to distinguish be- 
tween the various proposed gap symmetries. In previ- 
ous investigations for FeAs compounds the effect of non- 
magnetic impurities on the quasiparticle spectrum in the 
S stat e 26 ' 27 has been studied. Magnetic impurity ef- 
fects have so far only been discussed for the single band 
model with d x 2_ y 2 order parameter—, and for the two- 
band model for a classical local moment^ 8 .. At the same 
time, it is known that the local density of state around 
the magnetic impurity can provide significant informa- 
tion on the local electronic structure in the unconven- 
tional superconductor—. Note also that the influence of 



nonmagnetic 



,30,31 



and magnetic 32 impurities on the re- 



duction of superconducting transition temperature has 
been recently analyzed. 
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FIG. 1: (color online) Fermi surfaces of the two-band model 
(red and green thick curves) with different symmetry of the 
superconducting order parameters (nodal lines represented by 
the dashed lines), (a) extended s-wave symmetry, starting 
from fully gapped A s ±(k x ,k y ) (or A s ±(k x ,k y ) with a — 
0). With increase of the higher harmonics, a — 1,3,6,7,8 
in A s ± {k x , k y ) this gap becomes more anisotropic and fi- 
nally has accidental nodes on the electronic pocket around 
the M-point; (b) refers to the other extended nodal s-wave 
gap symmetry A s ±(k x ,k y ) with two separate nodal lines: 

A s ± (k x , k y ) for a = 1.2; a' = 0.15 and a = 1.17; a' = 0.08); 



(c)-(d) show d x 2_ y 2 
enforced nodes. 



and d xy gap symmetries with symmetry 



In this paper we investigate the effect of a single 
magnetic impurity on the local quasiparticle excitations 
around the impurity site. We use a minimal two-band 
model for the electronic structure of Fe pnictidcs which 
leads to the r(0, 0)- centered hole and M (it, 7r)-centered 
electron pockets. In section HT1 the Anderson model with 
a strong Hubbard repulsion for the localized f-electron at 
the impurity site, and a hybridization between conduc- 
tion bands and localized state will be introduced. We will 
treat this model in the infinite U limit where a slave boson 
representation may be used similar to Ref. [l^ where the 
Anderson impurity in the single band model with d x 2_ y 2 
order parameter has been studied. We then calculate 
the local density of states (LDOS) and discuss the signa- 
tures of possible Fe pnictide order parameter symmetry 
in its spectral and spatial characteristics. This quantity 
is accessible in STM tunneling spectroscopy^. The nu- 
merical results for the various cases will be discussed in 
section IIIII Finally in section IIVI we give a summary of 
our results and a conclusion. 
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FIG. 2: (color online) Calculated local density of states 
(LDOS) for the superconducting regime with various super- 
conducting order parameters: a) S ± : A s ±(k x ,k y ) and (b)- 



(f) Sf 



:(k x ,k y ) (with a = 1,3,6,7,8 respectively); for 



t t = -A /3, Vi = Vi = (green), Vi = V 2 = 0.5A (red), 
and Vi = V% = A (blue). 



II. THEORETICAL MODEL 

According to the band structure calculations^ as well 
as numerous ARPES results^ the Fermi surface topology 
of iron-based superconductors consists of the small size 
circular hole and elliptic electron Fermi surface pockets 
centered around the T— and (ir, 7r)-points of the folded 
BZ, respectively. The pockets are nearly of the same size 
which results in the nesting properties of the electron and 
hole bands at the antiferromagnetic wave vector, Q^p, 



i.e. £{j. = — £k + Q AF - Despite the fact that there are two 
electron and two hole pockets it has been argued^ that it 
is enough to consider only two of them (one electron and 
one hole pocket) because the instabilities of the two-band 
model are the same as in the four-band model. Following 
this suggestion we consider two bands that are given by 
diagonalized tight binding expression including hoppings 
up to the next nearest neighbors: 
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FIG. 3: (color online) LDOS for the superconducting state 
with various order parameters: A s ±(k x ,k y ) for (a = 

1.2; a' = 0.15) and (a = 1.17; a' = 0.08); for e f = -A /3, 
Vi = V 2 = (green), Vi = V 2 = 0.5A (red), and Vi = V 2 = 
A (blue). 



where creates an electron with spin a in band 7 

(7 = 1,2 refer to the electron and hole band, respectively) 
with wave vector k = (k x , k y ). The dispersion is then 
given in the tight-binding form^ 

tki = — 0.18 + 0.16 (cos k x + cos k y ) — 0.052 cos k x cos k y 



k k. 

£k2 = 0.68 + 0.38 (cos k x + cos k y ) — 0.8 cos cos 



(2) 



Here, £ki dispersion yields the hole Fermi surface pocket 
around the T point and £k2 gives the elliptic electron 
Fermi surface pocket around the M = (tt,tt) point, sec 
FigfT] The parameters have been chosen from the avail- 
able fit to the ARPES data^ (all in eV) and correspond 
to the hole doping of about 10%. The fl operator create 
the localized electron at the impurity site at the origin 
and U is its on-site Coulomb repulsion. Finally e/ is the 
f-band position, Vk 7 is the hybridization energy between 
localized electron and the conduction bands, and Ak is 
the singlet superconducting gap function. We choose val- 
ues for e/, and V kl (FigJ^J) such that the f-orbital is al- 
most filled (n,f = 1). 

Our model assumes the limit U — » 00 where doubly 
occupied f-states are excluded. This limit may be rep- 
resented by introducing the auxiliary boson b, with the 
constraint Q = b^b + flf a = 1— ■ In the mean field 
approximation (b = (b) = (b^)), the total Hamiltonian 
including the constraint is given by TCmf + A(6 2 — 1). 
Here TCmf is obtained as 



Hmf = + ^ k /M k )^k + (^ k / 3 2<?+ h.c), (3) 



where A is the Lagrange multiplier for enforcing the con- 
straint. The Nambu spinors are denoted by 

^1 = ( c kiT' c - k U' c k2T' c - k2 l)' and likewise $ = 



/i) f-f j /l)i while the matrices /3, are defined as 



A) = e/?o ® cr z ; 

A(k) 



1-C3 



(£k2CT 2 + A k a x ) 



/3 2 (k) 



l + ?3 



kr 



l-?3 



(4) 



Here Ui are the Pauli matrices acting in spin space, ?j 
are the Pauli matrices in the orbital space, and q <g> cjj 
denotes a direct product of the matrices operating on the 
4-dimensional Nambu space. Furthermore e.f = (e/ + 
A)/2, and Vk 7 = 6Vk 7 are effective hybridization and 
energy position of the impurity /-level, respectively. 

The local density of states (LDOS) near the magnetic 
impurity is obtained from analytic continuation iuj n — > 
E + i0 + according to 



N C (E, r) 



1 



Im [G c n (r,r,uj n ) + G§ 3 (r, r, w„)] ,(5) 



where u n = TrT(2n + 1) is the Matsubara frequency, and 
G c (r, r',ui n ) is a Fourier transformation of the matrix of 
the conduction electrons Green's function. 

The matrix Green's functions are defined as the 
imaginary-time ordered average 



G c (k,k';r) = 

G /c (k;r) = 

G c/ (k,r) = 

Gf(r) - 

Where G(r) = Tj2^G(Lo n )e 
regime TV (...) - 



(T r &(r)$,(0)>; 



(T r £(r)#(0)>; 

(T T 4(r)^ + (0)); 
(T r ^(r)^(0)). 



(6) 



At low temperature 



— » -F Im h~' d u. (...), 
where Z? 7 is bandwidth of conduction band 7. Using 
the standard equations of motion method, one can show 
that 



(iw n - /3i(k))G c (k, k', Wn ) = S ky + /3 2 (k)G /c (k',< 



(7) 



(iw„ - p )G fc (k,Lu n ) = ^/3 2 (k')G c (k',k, w „); (8) 



and 



(iw B -A,)G'(w n ) = l + ^/3 2 (k)G c /(k,^„); (9) 



(««;„ - /3i(k))G c/ (k, cj n ) = /3 2 (k)G / (w„). (10) 

Now using the equations above we find the full /-Green's 
function 



1 



■jj n 



(11) 



4 



where the /-self energy is given by 
£, = ]T/3 2 (k)Gg(k,^)/3 2 (k) 

k 



+(l->2; ft^-ft), 



(12) 



and the conduction electrons Green's function can be ob- 
tained byi2 



G c (k,k>„) = G c (k,u n )[5 k , k , +t(k,k';o;„)GS(k / ,a;„)]. 



(13) 



Here, Gg(k, w„) = (iw n — /3i(k)) - is the unperturbed 
Green's function of the conduction electrons, and the t- 
matrix is given by 



t(k,k';0 =/3 2 (k)G / ( W „)/3 2 (k') 



(14) 



In the following, the k dependence of the hybridization 
energy is neglected, i.e., we set Vk 7 = Vy, which yields 
that /? 2 (k) 

Minimization of the ground state energy with respect 
to b and the Lagrange multiplier A leads to the mean field 
equations 

A6 2 = E^4> 5>£ + 6 2 = l, (15) 

where the expectation values are defined by W(l a = 
(4ck T0 -), and n{ = (/]/<,). Therefore from Eq. {T 
we show easily that 

Afe2 = r l™o£ [ fl ( G n( k ' r ) ~ G 2 2 (k,r) 

k 

+V 2 (G{ 3 c (k,r)-G 4 / 4 c (k,r 



(16) 



and 



b = — lim 



G f ix(r) - G{ 2 {t) +G 3 / 3 (r) - G{ 4 (t)1 .(17) 



By solving the set of equations (TTTjl - lfrT)) . numerically one 
can find the values of if and b which are used as an input 
for the t- matrix.. 



III. NUMERICAL RESULTS 

We now discuss the results of numerical calculations 
for the central quantity N c (E,r) based on the previous 
analysis (obtained from Eqs. [TTIfl~7)l . In this section we 
focus on the energy and also spatial dependence of the lo- 
cal density of states (LDOS) (Eq. O at T = for various 
gap symmetries. As the main candidates we include dif- 
ferent types of the extended s-wave gaps which are fully 




FIG. 4: (color online) LDOS for the superconducting regime 
with order parameters: a) d x 2_ y 2 and b) d xy ; for e/ = — Ao/3, 
Vi = V 2 = (green), Vi = V 2 = 0.5A (red), and Vi = V 2 = 
A (blue). 



gapped on the hole pocket but possibly have accidental 
nodes on the electron pockets due to higher harmonics 
contributions. For completeness, we also include d-wave 
gap functions which have nodes on both electron and 
hole pockets. The latter however are not supported by 
ARPES result o 15 ' 38 ' 39 which suggest nodclcss gaps on the 
hole pockets. 

We begin our discussion with the anisotropic nodeless 
extended s-wave pairing function (S^). The Fermi sur- 
face illustration with respective node positions indicated 
by dashed lines is shown in FigfTJ(a). This gap function 
is given by 



A s ±(k x ,k y ) = — (cosfc^ 



cos ky ) , 



(18) 



and we show the result for the LDOS at the origin 
{r x ,r y ) = (0,0) in FigOJa) for A = 6 meV and e/ = 
—Ao/3, with different hybridization energies V\ = Vi = 
(green), Vi = V 2 = 0.5A (red), and V x = V 2 = A 
(blue). While the overall structure of the spectrum stays 
the same the increased hybridization leads to more pro- 
nounced bound state peaks within the gap. We observe 
that an increase of the hybridization energy, Vy, causes 
the position of peaks to move to higher absolute values 
of E, and a corresponding increase of their line width. 
By increasing the absolute value of the impurity energy 
level, e/ , the bound state also moves to higher absolute 
values of E. Notice that by restricting to the first order 
perturbation theory in effective hybridization (Born ap- 
proximation), we did not find a dramatic change in our 
results, for the chosen values of Vy. 

Furthermore, due to the different size of electron and 
hole pockets the onset of the continuum around \E/A\ is 
split into a double peak structure. For E > the lower 
and upper peaks correspond to hole and electron pockets, 
respectively. 

The addition of higher harmonics in the S x gap func- 
tion allows to tune the modulus of the gaps on hole 
and electron pockets independently while keeping the ba- 
sic property of having an opposite signs. As seen from 
Figda) this is equivalent to shift of the nodal line posi- 
tion closer to the electron pockets and increasing higher 
harmonic amplitude continuously one eventually produce 
an accidental node on this sheet. This gap function is 



5 



given by 



A s ± (k x , k y ) 



7- (cos k x + cos k v ) + a cos -A cos ^ 
2 2 2 



(19) 



and the result for the LDOS is plotted in Fig.(|2b)-(f) 
for different a parameters. At the symmetry points one 



has A 

' \o,o)| 



(0,0) 



A (l 



and A 



(ir,7r)) 



= -A . Then 



= a which also means that the dif- 



ference between the absolute values of the gaps on T- 
centered hole pockets and M-centered electron pockets 
increases with a. This can be clearly seen in Fig. HJb) 
(a = 1) where the peak at larger \E/Aq\ originating from 
the hole pocket is pushed to larger energies and for a = 3 
is no longer visible on the scale of Fig. [UJc). On the 
other hand the peak due to gap maximum on the electron 
pocket stays fixed around \E/Aq\ ~ 1 while a grows. At 
the same time a deep minimum and finally an accidental 
node of the gap develops on the electron pocket and leads 
to the increase of the low energy LDOS in Figs[2Jb)-(f). 
This type of low energy DOS may explain power laws for 
NMR relaxation rate and penetration depth observed in 
some pnictide s 20 ' 21 i 22 . On the other hand the position 
of bound state peaks caused by the magnetic impurity 
is apparently insensitive to the variation of a and the 
change of the underlying quasiparticle spectrum. 

However, the Sf gap function is not unique and ac- 
cidental nodes on the electron pocket may be obtained 
with a different type of modification with higher harmon- 
ics. This leads us to extended S^-wavc pairing described 
by another gap function 

A s ± (k x , k y ) = Ao [0.5(cos k x + cos k y ) 

+a cos k x cos k y + a' cos 4k x cos Ak y ] . (20) 

Its nodal structure is shown in (Fig. [2(b). Now there 
are two nodal lines, one located between the pockets 
which leads to an anisotropic but fully gapped order 
parameter on the hole pockets. The other accidental 
nodal line centered around the M point cuts the elec- 
tron pocket and leads to a finite low energy quasiparticle 
DOS as seen in the results of the Figs.([3](a) and[3jb) for 
a = 1.2; a' = 0.15 and a = 1.17; a' = 0.08 respectively 
The bound states due to impurity scattering appear again 
as pairs at similar energies as for the Sf case. 

Finally for completeness we also consider two sim- 
ple anisotropic d-wave order parameters, namely d x 2_ y 2 
(Fig. Q](d)) an d d xy (Fig. QJe)) which have symmetry 
enforced gap nodes. This leads to sign change of the 
gap function on the same FS pocket rather than between 
them. We consider the two candidates 



Ad x 2_ H 2 (^X) ky) 

Ad^y (k x , k y ) 



Ar 



(cos k x — cos k y ) , 



Aq sin k x sin k y . 



(21) 



In each case the nodal lines cross both Fermi surface 
pockets in contrast to the extended s-wave model. Wc 




o 
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FIG. 5: (color online) Density plot of the spatial distribution 
of the LDOS for the superconducting regime with extended 
s-wave order parameter, A s ±(k x ,k y ), for e/ = — Ao/3, Vi = 
V2 = Ao, and a) u> = ui r , b) ui = — u) r - (Red color correspond 
to large LDOS). 



note that current ARPES experiments have shown a fully 
gapped hole pocke t 15 i 38 i 39 though no experiments yet are 
available for P-based systems. 

The background of the LDOS is given by the typical V- 
shape of a d-wave order parameter. On top of it a single 
bound state peak due to the impurity scattering appears 
below the Fermi level. This is distinctly different from 
the extended s-wave case where always two bound states 
below and above the Fermi level appear symmetrically. 
A partly similar observation for a d x 2_ y 2 order param- 
eter with only a sing le sheet FS intended for cuprates 
was made in Ref. n» There one bound state peak was 
found for (r x ,r y ) along the anti-nodal direction and two 
peaks for the nodal direction. In our present d-wave case 
considered for the two-sheet FS of Fe pnictides the single 
peak appears for both nodal and anti-nodal directions. 

For clarifying of the position dependence of resonance 
peaks, Figj5] displays the spatial variation of the LDOS 
N c (u>,r x ,r y ) around the magnetic impurity at the reso- 
nance energies, (a) u = ui r and (b) u> = — ui r , for S 1 * 1 
gap symmetry with Ao = 6meV; e/ = —Ao/3, and 
V\ = V2 — Aq. It shows that the maximum amplitude of 
the LDOS appears close to the impurity site and decays 
non-monotonically further away from the impurity site. 
While the LDOS at u> r is rather isotropic the peak for 
— u> r shows a significant anisotropic LDOS in the plane. 
This anisotropy does not seem to result from special FeAs 
Fermi surface feature since it is also observed in the single 
parabolic band case in Ref. [H. 



IV. SUMMARY 

We have investigated the effect of magnetic impurity 
scattering in the FeAs pnictide superconductors. We 
used a simple two band model Fermi surface and cal- 
culated LDOS spectral and spatial dependence close to 
the impurity site for two types of extended s-wave super- 
conducting order parameters with inter-band sign change 
and for two d-wave order parameters with intra-band sign 



6 



change. In the former two impurity bound states ap- 
pear symmetrically around the Fermi energy at positions 
±uj r . The modulus of the bound state energy increases 
with hybridization strength V and impurity orbital en- 
ergy e/ monotonically. In the latter case only the lower 
bound state pole at — uj r appears in the LDOS for any 
direction from the impurity site. The background varia- 
tion of the LDOS is determined by the characteristics of 
the superconducting gap on the two FS sheets. The ex- 
tended S-wave order parameters may be tuned such that 
fully gapped behavior on the central hole sheet and acci- 



dental node structure on the zone boundary hole sheets 
appear naturally. In this case the spatial dependence of 
the LDOS for the two bound state peaks shows signifi- 
cant differences in the degree of spatial anisotropy. We 
conclude that the observation of two bound state peaks 
in tunneling experiments would be an important support 
for the extended s-wave gap function with intcrband sign 
change. The fine structure of the background continuum 
LDOS may give more detailed information on the type 
of the accidental nodal structure. 
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